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Abstract. In the present work, we investigated the correlation-induced localization-delocalization transition 
in the one-dimensional tight-binding model with fractal disorder. We obtained a phase transition diagram 
from localized to extended states based on the normalized localization length by controlling the correlation 
and the disorder strength of the potential. In addition, the transition of the diffusive property of wavepacket 
dynamics is shown around the critical point. 

PACS. 72.15.Rn Localization effects - 72.20.Ee Mobility edges - 71.70.-|-h Metal-insulator transitions - 
71.23.An Theories and models;localized states 


1 Introduction 

In some kinds of one-dimensional Schrddinger operators 
with random potential, the eigenstates are mathematically 
proven to be localized wim- The one-dimensional tight- 
binding models with an ergodic and stationary random 
potential have positive Lyapunov exponent of the wave- 
function with probability 1 (G-M-P theorem) [3]. The ex¬ 
istence of the positive Lyapunov exponent is necessary and 
sufficient condition for a pure point set spectrum of the op¬ 
erators, and then all the eigenfunctions exhhibit the expo¬ 
nentially decay in the thermodynamic limit. Kotani’s the¬ 
ory states that if the potential sequence is nondeterminis- 


Shima et al and Kaya also showed that the IDDS with 
FFM potential have a critical disorder strength Wc sepa¬ 
rating the conducting and insulating phases, and the Wc is 
independent of the spectrum index a(> 2) [1^113) . Peter- 
cen and Sandler insisted that effect of the anticorrelation 
is also important to understand the transition due to the 
correlation of the sequence [23] . 

On the other hand, the LDT due to the differentiabil¬ 
ity of the potential function also exists without contradic¬ 
tion with the Kotani’s theory. Very recently, Garcia and 
Cuevas studied the transition based on the differentiability 
of the disorder potential as a necessary condition for the 
delocalization |20l|2T]. A certain degree of differentiability 


tic under the following conditions, (i)stationarity, (ii)ergodici^sures that the potential in neighboring sites is strongly 


(hi) integrability, then there is no absolutely continuous 
(a.c.) spectrum of the operators [5]. These theorems can be 
proven true for continuous and discrete one-dimensional 
disordered systems (IDDS) |6l|7l|8]. However, the neces¬ 
sary and sufficient condition for the exponential localiza¬ 
tion has not been found yet. 

There is a possibility that the correlation effect of the 
potential sequence breaks the strong exponential local¬ 
ization and generates a localization-delocalization tran¬ 
sition (LDT) in the IDDS. Indeed, many authors numer¬ 
ically observed the correlation-induced LDT by using the 
some potential sequences with power spectrum S{f) ^ 
I//“(a > 2) as a potential by Fourier filtering method 
(FFM), where / denotes frequency [S HTUlITTlIT^IIdllldlllS] . 
Such a potential obviously breaks a necessary condition 
(i), i.e. it is nonstationary. Note that the non-stationarity 
can be satisfied for a > I. In addition, the LDT with 
mobility edges, has been found numerically for the one¬ 
dimensional tight-binding model with the FFM model as 
the potential strength W decreases [TBl[r7l[T51ll91l20l[^ 

na. 


correlated. They modeled the sequences with power-law 
spectrum by Weierstrass function with fractal dimension 
D. As a. result they also numerically suggested that the 
transition takes place at the critical value Dc = 3/2 by 
means of the distribution of the energy level-spacing in 
the weak disorder limit. However, studying the effect of 
the disorder strength in the IDDS with the Weierstrass- 
type potential has been still unclarified. 

In this paper, we have numerically studied the correlation- 
induced localization-delocalization transition by using tight- 
binding model with Weierstrass potential used by Garcia 
and Guevas m- The finite-size scaling analysis for the 
normalized localization length at band center numerically 
suggests the existence of the transition around D ~ 3/2 
independent of the potential strength in the relatively 
weak disorder regime. On the other hand, in the relatively 
strong disorder regime, the critical fractal dimension Dc 
becomes smaller value than D = 3/2 dependently on the 
potential strength. Furthermore, we have investigated the 
quantum diffusion of the initially localized wavepacket in 
the system. The transition from the localized state to bal- 
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listic states occurs around I? ~ 3/2 without scale invariant 
subdiffusive behavior. 

This paper is organized as follows. In the next section, 
we introduce the IDDS with the Weierstrass potential and 
some preliminary calculations. In SectOl we present global 
behavior of the IT—dependence and II—dependence in the 
LDT by the numerical calculation of the normalized lo¬ 
calization length. A phase diagram in the D — W space is 
also given. In SectUl we show that the dynamical prop¬ 
erty of the system ought to be represented by the time 
dependence of the degree to which the initially localized 
wavepacket would be spread. Summary and discussion are 
presented in the last section. In an appendix, we give the 
behavior of the autocorrelation functions for the potential 
sequence. 


2 Model and preliminary calculation 

2.1 model 


We consider the one-dimensional tight-binding model de¬ 
scribing single-particle electronic states in the site repre¬ 
sentation as 


u{n — 1) + u{n + 1) + WV {n)u{n) = Eu{n), (1) 

where E and {u{n)}^^Q are the energy and state of the 
system, respectively. The {V{n)}^^Q and W are the on¬ 
site energy sequence and the strength, respectively. To 
model the correlated and non-differential disorder poten¬ 
tial for V{n){n < N) in Eq.(IT|), we use the following form: 


L 


E(n) = C'^ 


sm{27ra^n/N + (pk) 

g^{2 — D)k 


( 2 ) 


where a is a constant value (a > 1) related the scale- 
invariance and D is a fractal dimension (1 < D < 2). 
{'/fe}fc=o 9-^® random independent variables chosen in the 
interval [0,27r]. C is the normalization constant which is 
determined by a condition 

\/< V{nY >-< V{n) >2 = 1, (3) 

where (...) indicates the average over realization of the 
phases in Eq.(I2|). 

If we set n/N = x, (j)k = 0, the potential sequence be¬ 
comes ’’Weierstrass function” with continuous and indif- 
ferentiable everywhere by taking a continuous limit N —^ 
oo and L oo. Therefore, the potential will be shortly 
transfered to as ’’Weierstrass potential” in this paper, and 
we set a = 2 and L = 50 thorough this paper without loss 
of the generality and accuracy of the numerical calcula¬ 
tion. The power spectrum S{f) of the Weierstrass func¬ 
tion is empirically characterized by the fractal dimension 
D as, 


S{f)- 


1 


( 4 ) 


In compassion with the form S'(/) ^ I//“, 


D = l + 


3 — 0 
2 


( 5 ) 


Note that the condition a > 2 for the LDT in the FFM 
potential lilllllS] corresponds to a condition D < 3/2. 
Increasing a corresponds to increasing correlations up to 
long-range correlated disorder. The analytical property of 
the autocorrelation function for the Weierstrass potential 
V (n) is given in appendix It is found that the correla¬ 
tion linearly decreases from 1 for D = 3/2. 

In addition to the long-rang correlation, the fractal di¬ 
mension D also controls the degree of the differentiability 
of the potential function. The degree of the differentia¬ 
bility increases along with the decrease of the fractal di¬ 
mension D. The smoothness of the potential fluctuation 
can also induce the delocalization of the quantum states, 
which property is directly related to analyticity of the po¬ 
tential function in the continuum limit, as pointed out by 
Garcia and Cuevas EOIIII]. They have numerically found 
that LDT at Dc = 3/2 for the sufficiently weak disorder 
regime by using the nearest-neighbor level-space distribu¬ 
tion of the energy spectrum. The result is not at odds with 
Kotani’s theory because the Weierstrass potential become 
non-stationary for 1 < D < 3/2. 


2.2 preliminary calculation 


The finite size Lyapunov exponent » 1) of the 

one-dimensional systems can be defined by 


In = 


ln(|u(iV)|2 + |u(iV+l)|2 
2N 


( 6 ) 


with initial state m( 0) = u(l) = 1. Then the Lyapunov 
exponent 7 and the localization length ^ are given by 
7 = lim„_>.oo 7Ar and by ^ = I/ 7 , respectively. The en¬ 
ergy dependence of 7 is strongly correlated with the den¬ 
sity of states that can be obtained from some experiments 
for real materials. In addition, we define the normalized 
localization length (NLL), 


_ m) 1 

N < 7 Ar > W’ 


( 7 ) 


, where (...) denotes the ensemble average. It is useful to 
study the LDT that Apf decreases (increases) with the sys¬ 
tem size N for localized (extended) states, and it becomes 
constant for the critical states. 

Figure (IDa) shows the energy dependence of the Lya¬ 
punov exponent 77 V for some values of the fractal dimen¬ 
sion D. The Lyapunov exponent 'Jn at the band cen¬ 
ter F = 0 decreases as the value of the D decreases. It 
suggests a possibility of the delocalized states (extended 
states) at the band center E = 0 for small values of 
D{< 3/2). Figure[TJb) and (c) show the potential strength 
dependence of the Lyapunov exponent < jn > and the 
renormalized localization length Ajv at the band center 


p-2D- 
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E = Q for some parameter sets. In the weak disorder 
regime W «1, the numerical data lead to 

In oc W, ( 8 ) 

despite of the system size, which is equivalent to case of 
the uncorrelated disordered system. In the strong disorder 
regime, the Lyapunov exponent does not depend on the 
system size N. Figure (TJc) shows the IF—dependence of 
the NLL Ajv. It is found that the localization length is 
larger than system size in the week disorder limit. The re¬ 
sults suggest that for the strongly correlated limit, there 
is a possibility of the states with < yjv >—?> 0 and/or 
> 1 in the thermodynamic limit iV —> oo, which cor¬ 
respond to a delocalized phase with the extended states in 
week disorder limit. In the following section, we numeri¬ 
cally investigate the U—dependence, IF—dependence and 
IV—dependence of the NLL in the IDDS with the potential 
©• 



Fig. 1. (Color online) (a)The Lyapunov exponent yjv as a 
function of energy E for a weak disorder case, IF = 0.2 and 
varying the fractal dimension D. The system and ensemble 
sizes are N = 2^^ and 2^°, respectively. (b)The Lyapunov ex¬ 
ponent 7 ]v and (c)normalized localization length An at the 
band center F = 0 as a function of the potential strength IF 
for the several values of the D. This system and ensemble sizes 
are N = 2^'^ and 2^'^, respectively. 


3 Localization-Delocalization Transition 

This is the main section of the present paper. In what 
follows, we investigate the NLL at band center E = Ohy 
changing the system size for some typical parameter sets 


{W,D). The typical basis size N and ensemble size used 
here are TV = 2® ~ 2^^ and 2^° ~ 2^^, respectively. The 
robustness of the numerical calculations with respect to 
the system size has been confirmed in each case. 


3.1 H-dependence and IF-dependence in the wide 
parameters space 



Fig. 2. (Color online) The normalized localization length An 
as a function of the potential strength IF with a fixed value of 
D for N = 2l^2l^2l^2l'‘,2l^2l®. (a)D = 1.7, (b)D = 1.5, 
(c)D = 1.4, {d)D = 1.3. Note that the data are plotted in 
double-logarithmic scale. 


Figure [5] shows IF-dependence of the NLL for some 
values of the D in the relatively wide range. The IF-dependence 
relatively smoothly drops down around IF ^ 1 in the 
same way for all cases with different system size. There 
exists the strong system size dependence for the strong 
disorder regime IF >1. Furthermore it is found that 
there is a quite different feature between the cases for 
D < 3/2 and ones for D > 3/2. In the case oi D = 1.7, 
the curve of the IF-dependence for the different system 
size (TV = 2^^ ^ 2^^) do not cross within this regime 
of IF with each others. This fact implies that all states 
around E = 0 are localized in the thermodynamic limit 
and there is no transition for changing the value of IF. 

On the other hand, in the cases of D < 3/2, the system 
size dependence for weak disorder regime IF < 1 becomes 
very weak, and the IF-dependence sharply decreases at 
certain value of IF. Apparently, we can expect that for 
D < 3/2 the IF-dependence of the NLL shows a sharp 
jump in the thermodynamic limit TV —>■ oo. This feature 
suggests the existence of a transition to delocalized states 
in the limit TV —>■ oo. 
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Fig. 3. (Color online) The normalized localization length 
yliv as a function of D with a hxed value of W for N = 
2 ii 212 213^ 214^ 215, 2i®. (a)VF = 0.9, (b)fF = 0.8, (c)W = 0.7, 
{d)W = 0.6, {e)W — 0.4, {i)W = 0.2. Note that the vertical 
axis are in logarithmic scale. 


In Figl3] we show the dependence of the NLL for 
some values of W in order to find out the critical value 
Dc that all curves of different system size intersect at the 
value. At least, the existence of the intersection corre¬ 
sponding to a transition point can be observed around 
D = 3/2 for the relatively weak disorder strength W < 
0.7. 

In the following two subsections, we investigate the 
details of the LDT with focusing on the relatively weak 
disorder regime {W < 0.7) and the relatively strong dis¬ 
order regime (0.7 < W < 0.8), respectively. 


3.2 Weak disorder regime 

Figure SKa) shows the system size dependence of the NLL 
for a relatively weak disorder strength W = 0.4. It is found 
that the Af—dependence changes the decreasing function 
to the increasing function as the fractal dimension de¬ 
creases, and the A^—dependences are algebraic. Generally, 


the quantum states can be classified by the exponent S of 
the Af—dependence of the NLL Ajv when it behaves as, 

Ajv - Af^ (9) 

The exponent, 5 < 0 for the localized states, 5 > 0 for the 
extended states, and 5 = 0 for the critical states. The ex¬ 
ponents 6 obtained by the least-square method are shown 
in Figl^ The exponent 6 decreases with respect to the 
fractal dimension D as 

5 ~ (3/2-D)-3 ®®. (10) 

It seems that in the Af—dependence of the NLL the sign 
of the index 5 changes from negative to positive one at the 
point Dc = 3/2 for W = 0.4 in FiglJKa). 

As a result it is suggested that the LDT takes place 
around the transition point Dc = 3/2 independent of the 
disorder strength in relatively weak disorder regime W < 
0.7, as shown in FiglU Figl^l and FigH] 

However, it should be noticed that the transition does 
not obey the standard one-parameter scaling theory (OPST) 
of the localization [M] if the A^—dependence seen in Fig|3Ka) 
continues in the thermodynamic limit N —> oo. In the 
present stage, it is difficut to get the A^—dependence with 
adequate accuracy for the larger system size because of 
the limitation of the numerical calculation. 


3.3 Strong disorder regime 

We examine the Ajy behavior to get the delocalization in 
the relatively strong disorder regime {W > 0.7). 

Figure HKb) shows the system size dependence of the 
NLL for W = 0.98. It clearly shows that the wavefunc- 
tion goes to localized states in the thermodynamic limit 
as Atv ^ A^“3 irrespective of the fractal dimension for 
the relatively large disorder strength W = 0.98. Surely we 
have to investigate the smaller values of D{^ 1) to make 
the delocalized states for W = 0.98 although it is very 
hard. 

Figure [5] illustrates the more detailed behavior of 
as a function of D and W for D <?> 12 and the relatively 
strong disorder W > 0.7. It is found that the crossing 
points gradually shift and converge to certain value in the 
thermodynamic limit N ^ oo. This result suggests that 
the critical value of D decreases less than 3/2 for the rel¬ 
atively large potential strength W. 

The A^—dependence of the NLL for FiglOKa) and (d) 
are drawn in FiglTKa) and (c). The latter both reveal a 
clear change from the decreasing function to increasing 
one, as the D and W decreases, respectively. Next, we try 
to construct the scaling function by the parallel shift of 
the horizontal axis. If the OPST around the LDT exists, 
then the NLL behaviors as 

An = f{N/0. ( 11 ) 

where /(...) is the scaling function and ^ is the amount of 
the parallel shift, corresponding to the localization length 
or correlation length. In FiglTKb) and (d), we show the 
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Fig. 4. (Color online) The normalized localization length yljv 
as a function of the system size N with several values of the 
fractal dimension for the potential strength, {a)W = 0.4 and 
{h)W = 0.98. Note that the data plotted in double-logarithmic 
scale. 



Fig. 5. (Color online) The index 5 of Eq.® as a function 
of the fractal dimension D for the case W = 0.4 in FigU 
The result has been obtained by the least square fit for the 
Y—dependence of the normalized localization length ylAr. The 
straight line with a slope —1.88 is shown as a guide for eye. 


result for the latter shift in the FiglTKa) and (c), respec¬ 
tively. The general forms of the graph are similar to the 
results for IDDS with FFM-potential [TH]. 

The upper (lower) brunches are delocalized (localized) 
regimes. As seen in FiglT^b) and (d), it seems that the 
data of the lower branches look like being on a common 
curve while the data for the upper branches do not belong 
to an universal curve, as far as present data show. The 
A^—dependence oi with adequate numerical accuracy 

for the larger system size is necessary to perform the finite- 
size scaling analysis around the critical point. 



Fig. 6. (Color online) Some enlargements of the 
D—dependence and VF—dependence of the normalized 
localization length yljv. (a)VF = 0.75 and {h)W = 0.8, 
{c)D = 1.2, {d)D — 1.3. The system size N is varied as 2® 
~ 2^^. The averaging is taken over 2^^ realizations. Note that 
the data are plotted in the semi-logarithmic scale. 


3.4 phase diagram 

We determined the critical values of D and/or W by the 
U—dependence and W—dependence of NLL as seen in 
FigI21 FigOl 

FigH] andFigini Figure [5] shows the phase diagram in 
the D — W space separating the localized and delocal¬ 
ized states. The most interesting point here is that the 
critical disorder strength Wc depends on D in the strong 
disorder regime {W > 0.7). This result implies that the 
LDT with the H—dependent critical disorder strength is 
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Fig. 7. (Color online) The normalized localization length An 
as a function of the system size N with (a)W — 0.75 for the 
different values of Z). /Ijv as a function of the system size N 
with {c)D = 1.2 for the values of W. Panels (b) and (d) are 
the scaling function construction from the numerical data in 
panels (a) and (c), respectively. Note that all are plotted in 
double-logarithmic scales. 



Fig. 8. (Color online) Phase diagram of the localization- 
delocalization transition in the W — D plane. The transition 
points (denoted by the filled blue squares) have been obtained 
by numerical plots of the normalized localization length as 
shown in FigO Fig[3]and Fig[6l The inset shows the diagram 
in the D — W plane. The open symbols (circles, triangles, rhom¬ 
buses) indicate the parameter sets used in the quantum diffu¬ 
sion in Fig[5]and FigfTIJl 


different from the one observed in the IDDS with FFM- 
potential, in which it is Wc = 2 in our notation indepen¬ 
dent of the spectrum exponent a [18) . And our result in 
part coincides with one for the FFM-system by Kaya that 
the critical values of W increases when a increases, while 
the critical exponent v decreases when a increases |13j . 
Accordingly this result suggests that the phase diagram 
in the a — W space might be different even if the poten¬ 
tial sequences are characterized by the same exponent a 
of the power spectrum in the LDT in the IDDS. 


4 Quantum diffusion 

It can be expected that for sufhciently differentiable po¬ 
tentials a band of delocalized states occurs due to destruc¬ 
tion of the interference effects in the reflected components 
of the wavepacket. 

In this section, we examine the quantum diffusion of 
the initially localized wave packet by changing the param¬ 
eter D. We monitor the mean square displacement (MSD) 
of the wavepacket, 

™2(f) = no)^\u{n, f)p, (12) 


where uq is an initially localized site. The quantum time- 
evolution is given by, 

ih — ^ ^ = u{n -I-1, t) -|- u(n — 1, t) -|- WV{n)u{n, tjj)-3) 

at 

where n = 1,2,..., N with initial state u(n, f = 0) = 5n,no 
and = 1. 

Figure [S] shows the time-dependence of the MSD for a 
fixed potential strength W = 0.8 and varying the fractal 
dimension D. The parameter sets we used are denoted in 
Figl8]by some open symbols. Apparently, the m2 ballisti- 
cally grows for the relatively small values of D, and it is 
localized for the cases with D > 3/2. 

Next, FiglTO] shows the diffusive properties for some 
values of the potential strength with fractal dimension, 
D = 1.4, D = 1.6, denoted by the open symbols in FiglH] 
It follows that the ballistic-like motion {m2 , {3 > 1) 

can be obtained for all cases with D = 1.4, while it is well- 
localized for the relatively larger value of W in the cases 
with D = 1.6. The all of the cases oi D = 1.6 are localized 
for the long-time calculation. (It is not shown here.) 

As a result, it seems that in the quantum diffusion the 
critical value Dc also around D ^ 3/2 which is consistent 
with the result by the normalized localization length An in 
the last section. However, clear subdiffusive behavior have 
not been observed around the critical point Dc = 3/2, 
which is different from the localization-delocalization tran- 
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Fig. 9. (Color online) The second moment m 2 as a function 
of time for some fractal dimensions with W = 0.8. The bold 
curve corresponds to the critical case D = 3/2. We set h = 1 
and St = 0.05. The system size and sample sizes are N = 2^® 
and 10, respectively. 




Fig. 10. (Color online) The second moment m 2 as a function 
of time for some values of W {W = 0.6, 0.7, 0.8, 0.9) at (a)D = 
1.6 and {h)D — 1.4. We set h = 1 and St = 0.05. The system 
size and sample sizes are N = 2^® and 10, respectively. The 
inset shows the log-log plots. 


sition in three-dimensional disordered systems and poly- 
chromatically perturbed IDDS, in which cases an asymp¬ 
totic behavior, m 2 ~ has been obtained at the critical 
point This discrepancy is not so surprising if 

the correlation-induced LDT deviates from the standard 
OPST as shown in the previous section. 


5 Summary and Discussion 

In this work we numerically investigated the combined 
effects of the disorder strength W and the order of the 
long-range correlation characterized by fractal dimension 
D in the one-dimensional tight-binding model with the 
Weierstrass potential. The result based on the normalized 
localization length strongly suggests that quantum states 
are localized for D > 3/2, whereas we have obtained crit¬ 
ical disorder strength Wc which separates extended and 
localized regimes for D < 3/2. In particular, the result 
that the critical value Dc = 3/2 is not depend on the W 
is consistent with the result given by Garcia and Cuevas 
in the weak disorder limit. While the critical value of the 
fractal dimension depends on the potential strength in 
the strong disorder regime. Furthermore, the localization- 
delocalization property reflected on the quantum diffusion 
of the initially delocalized wavepacket although the subd¬ 
iffusion could not be observed at the critical case. The dis¬ 
order strength dependence of the correlation-induced tran¬ 
sition is one of the interesting features in the localization- 
delocalization transition. 

Concerning the power spectrum S{f) ~ !//“,« = 2, 
the localization-delocalization transition point Dc = 3/2 
is consistent with one predicted in the system with some 
other correlated potential such as FFM potential IMM 
[Miisniisi]. Moreover, there are the other types of long- 
range correlated potential with discrete values such as bi¬ 
nary and ternary sequences [2911^ . One of the common 
points between the discrete and continuous models is non¬ 
stationary of the potential sequences caused by the long- 
range correlation. Accordingly, the indifferentiable every¬ 
where condition may be unihed into the nonstationary 
condition (a > 2) for the delocalization in the IDDS. 

An interesting question to ask is whether it is pos¬ 
sible to characterize the correlation-induced localization- 
delocalization transition by the standard one-parameter 
scaling theory. The study of the hnite-size scaling and 
critical exponents with adequate numerical accuracy are 
future challenge. 


A Correlation function 

The normalized autocorrelation function C(a, n, m) of the 
Weierstrass potential sequence V{n) can be analytically 
calculated, as given for the FFM-potential by Petersen 
and Sandler [55]. The explicit form becomes, 

< V(n)V{m) > 

< V(nY > 


C(a, n, m) 


(14) 
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L/2 


^„-2(2-Z3)fccos 




27ra 


u \n — m\ 


L/2 

/c=0 


(15) 


-2(2-D)k 


where < ... > denotes the ensemble average over the in¬ 
dependent phases {cjik} in Eq.©. We set the distance 
between positions £ = n — m, and impose the periodic 
boundary conditions on the correlation function defined 
in £ G [0, L/2]. In addition, we set r = 2£/L and r G [0,1] 
for the thermodynamic limit L —>■ oo. Then the autocor¬ 
relation function can be written as 


C{a,r) 


OO 

cos(7ra^r) 

k^O _ 

OO ’ 

a-2(2-r>)fe 

fc=0 



(16) 



In the critical case, D = 3/2, it becomes 

OO 

cos(7ra^r) 

C{a,r) = ^^^ -. (17) 

E«-'^' 

fc=0 

Figure [TT] shows the autocorrelation function C{a = 
2,r) given by Eg. (1161) for various values of some fractal 
dimensions. It follows that the correlation function rapidly 
decays with complex fluctuation for D = 1.9. Moreover, 
the correlation function becomes concave for H < 3/2 
and it linearly decreases near r ~ 0 for the critical value 
Dc = 3/2. The smaller the fractal dimension D becomes, 
the correlation becomes more negative. The inset of the 
Fig. [TT] shows the values of the correlation function at 
r = 1 as a function of the fractal dimension. It is noted 
that the correlation function goes negative value at the 
thermodynamic limit r = 1 for D < 3/2. Such a property 
has been pointed out by Petersen and Sandler in the case 
of the FFM-potential [23] . 
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